Ambiguities on the Hamiltonian formulation 
of the free faUing particle with quadratic 

dissipation 

G. Lopez^, P. Lopez^, X. E. Lopez^, and R. I. Castro^ 

^ Departamento de Fisica de la Universidad de Guadalajara 
Apartado Postal 4-137 
44410 Guadalajara, Jalisco, Mexico 

2 Facultad de Ciencias de la UNAM 
Apartado Postal 70-348, Coyoacan 04511 Mexico D.F. 

February 2, 2008 

Abstract 

For a free falling particle moving in a media which has quadratic velocity force 
effect on the particle, two equivalent constants of motion, with units of energy, two 
Lagrangians, and two Hamiltonians are deduced. These quantities describe the dynamics 
of the same classical system. However, their quantization and the associated statistical 
mechanics (for an ensemble of particles) describe two completely different quantum and 
statistical systems. This is shown at first order in the dissipative parameter. 

PACS: 03.20. +i, 03.30. +p, 03.65.-w 

1 Introduction 

It is well known that the Lagrangian (therefore the Hamiltonian) formulation for some systems 
of more than one dimension may not exists (Douglas 1941). Fortunately for our study of the 
nature up to now, most of our physical systems have avoided this problem, and the whole 



quantum and statistical mechanics of non-dissipative systems can be given in terms of a 
Lagrangian or Hamiltonian formulation. Now, for dissipative systems there have been two 
main approaches. The first one consists of keeping the same Hamiltonian formalism for the 
whole system where the interacting background is included, as a result, one brings about a 
master equation with the dissipation and diffusion parameters appearing as part of the solution 
(Caldeira and Legget 1983, Unruh and Zurek 1989, and Hu et al 1992). This approach has 
its own merit, but it will not be followed in this paper. We will follow the second approach 
which consists in to obtain a phenomenological velocity depending Hamiltonian, representing 
a classical dissipative system, and to proceed to make the usual quantization (or statistical 
mechanics) with this Hamiltonian. 

Within this last approach, one can, additionally, study the mathematical consistence of 
the the Hamiltonian formalism in quantum and statical mechanics. It is also known that even 
for one-dimensional systems, where the existence of their Lagrangian is guaranteed (Darboux 
1894), the Lagrangian and Hamiltonian formulations are not free from problems (Havas 1973, 
Okubo 1980, Dodonov et al 1981, Marmon et al 1985, Glauber et al 1984, "Lopez 1998, and 
''Lopez 1999). One of the main problems is the implication on the quantization of the associ- 
ated classical system when different Hamiltonians describe the same classical system (^Lopez 
2002). This ambiguity has already been studied for the harmonic oscillator with dissipation 
and some general system ('^Lopez 1996). In this paper, we want to show explicitly this am- 
biguity by studying the free falling particle within a medium which has the effect on the 
particle of producing the dissipation. This dissipation depends quadratically on the velocity 
of the particle. Firstly, two constants of motion are deduced for this system. Secondly, with 
these constants of motion two Lagrangian and two Hamiltonian are obtained using a known 
procedure ('^Lopez 1996 and ^Lopez and Hernandez 1989). Finally, using the Hamiltonian 
expression at first order in the dissipation parameter, the resulting eigenvalues of their as- 
sociated quantum Hamiltonian and their associated statistical mechanics properties (for an 
ensemble of particles) are shown. 
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2 Constants of Motion 



The motion of the particle of mass m falhng under a constant gravitational force, —mg, where 
g is the constant acceleration due to gravity, which is within a dissipative medium which has 
the effect on the particle of producing a force proportional to the square of the velocity of the 
particle, ax^ for i; < 0, can be described by the following autonomous dynamical system 

dx dv a ^ , 

-n=^^ -77 = -9 + -V : (1) 

at at m 

where the variable x represents the vertical position of the particle, and v represents its 
velocity. A constant of motion for this system is a function K — K{x, v) such that dK/dt — 0, 
i.e. it satisfies the following equation (•''Lopez 1999) 

OK ( « n\ dK , . 

The general solution of this equation is given by (John 1974) 

K^{x,v)^G{C{x,v)) , (3) 

where G is an arbitrary function of the characteristic curve C{x, v). This characteristic curve 
can be given in two different ways as 

+ 9^ , (4) 

za \ mg j 

or 
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Considering that one must obtain the usual constant of motion (Energy) expression for a 
equal to zero, the functionality of G in Eq. (3) is determined for each above characteristic 
(G(Ci) = mCi and G{C2) = —{mg/2a)C2 — m'^g/2a), and the following constants of motion 
are gotten 

v) = In (l - —vA + mgx , (6) 

2a \ mg ) 

and 

Note that the following limit is gotten 

lim K^^ = -mv'^ + mgx i = 1, 2 . (8) 

a— >0 2 



3 



3 Lagrangians and Hamiltonians 

Using the know expression (Kobussen 1979, Leuber 1987, Yan 1981, and '^Lopez 1996), 

L{x,v)=v J ^^^dv , (9) 
to get the Lagrangian, the Lagrangian associated to Eq. and Eq. ([7]) are 



and 



p)^' = m\j arctanh \ J v ) (12) 



if«(x,p) = -^ln 



1 - tanh^ 



a p 
mg m 



and 
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(10) 



L^^Hx, v) = m\l^^ V arctanh ( . /-^ f | H In ( 1 v'^ ] — max 

\ a \ \ mg J 2a \ mg ) 

and 

4^n-,^^) = ^(^ + -^^0e---/--^. (11) 
la \ m J la 

Their generahzed hnear momenta [p = dL/dv) are 



= ^^e-'""/™ . (13) 
Thus, their associated Hamihonians, H{x,p) = K{x,v{x,p)), are given by 



+ mgx (14) 



H!^\x,p) = !^(-g + ^e^"-/- + ^ , (15) 

la \ m-^ J la 

where one has made the substitution of pa"* and pa"* by just p. One must note that the 
following limits are gotten 

hm L^*-* (x, v) = —mv^ — mgx , (16) 

a— >0 2 

hmp^*^ = mv , (17) 
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and 

^2 



r 

In 

At first order in the dissipation parameter a, the Hamiltonians are given by 



hm Hj^\x,p) = h mgx i = 1,2 . 

a^o 2m 



and 



2 

H^^\x,p) = h mgx — (19) 

2m 12m^g 

H^^)(x,p) = — + mgx + - gxA . (20) 
2m \ m^ J 

4 Quantization at first order in perturbation theory 

Our Hamiltonians Eq. (|T9l) and Eq. fl20|) can be written as 

H^'\x,p) = Ho{x,p) + W^'\x,p) for i = 1,2 , (21) 

where Hq represents the non dissipative part of the Hamihonian, 

p2 

Ho{x,p) = - ^mgx, (22) 

2m 

and W^^^ represents the contribution of the dissipation at first order in a, 

W'^'\x,p) = --^^ (23) 
12gm^ 



and 



W^^\x,p) = a(^-gxA 
\ m^ J 



(24) 



The Schrodinger equation, ihd'^{x, t)/dt = H^'^\x,p)'if{x, t) represents an stationary problem. 
Therefore, in order to get the quantization of the system, one just need to solve the following 
eigenvalue problem 

i^«$„(a;) = ^«$„(a;) , (25) 
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where H^''^ is the associated Hermitian operator of Eq. fl21l) . Of course, one must not allow the 
particle to go beyond down the surface level. Thus, Eq. (1221) is representing the Hamiltonian 
of the quantum bouncer (x > 0) (Gean-Banacloche 1999), where the eigenvalue problem 

^o(£,p)^f (a^) = 4°Vi°n^) (26) 
has the eigenvectors and eigenvalues solution given by 

(0) Aijz - Zn) , . 

and 

= mglgz^ . (28) 

The functions Ai and Ai' are the Airy function and its differentiation respect to z. The 
variable z is defined as 2; = x//^, where Ig is given by Ig = (h? /2m?gY^^, and Zn is the nth-zero 
of the Airy function {Ai{—Zn) =0). In fact, the bouncing problem has already been studied 
for linear and quadratic dissipation (^Lopez 2004). For the later, the correction given to the 
eigenvalue problem using Eq. at first order in perturbation theory is 

{n\W^^^\n) = a-^^ . (29) 
15 

Now, using the relation {n\d^ / dz'^ln) = z'^/5, the correction at first order in perturbation due 
to Eq. (1251) is given by 

fc4 2 

/n\wW\n) = -a—^ . (30) 

Therefore, for the same classical dynamical system we have two different associated quantum 
systems which have completely different quantum dynamics, which is shown through the 
eigenvalues 

E^^ = + a^-^ (31) 



and 



= E(o) _ a-^^ . (32) 
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5 Classical Statistical model for dissipation 



Consider a system of iV = A'^i + particles, where A'^i particles are small of mass mi, and 
N2 particles are are big of mass m2 (m2 3> mi). The small particles move under the action 
of an external force with components (0, 0, —mg) and suffer collisions with the walls of the 
container which consists in a narrow-square shape pipe of cross sectional area L^. In addition, 
each small particle can have occasional (stochastic) collision with the big particles, when they 
are added, establishes the dissipative medium where the big particles will move. The big 
particles move in this dissipative medium, and it is assumed that, since this type of collision 
does not occur frequently, its average effect may have neglected contribution on the dynamical 
macroscopic variables of the system. Newton's equations of motion for this system can be 
written as 

rriiqiij =0 j = 1, . . . , Ni, i = x,y (33) 
miQuj = -mig j = 1, . . . , iVi (34) 
m2q2ik = a {q2ikf k = 1, . . . ,N2; i = x,y (35) 

m2q2zk = a {q2zkf - m2g k=l,...,N2, (36) 

where qaij, qaij and (jaij are the generalized coordinates, velocities and accelerations of the 
light-small (a = 1) and heavy-gross (a = 2) particles, and the parameter a characterizes the 
dissipative medium. The Hamiltonian associated the the motion of 1-particle, Eq. fl5^ and 
Eq. ([MD, is given by (''Lopez et al 1997) 



Hl-x,y,z = 

The Hamiltonian associated to Eq. (!35|) is written as 



3 „2 

EPiij , 

1=1 



(37) 



^222 /O \ 

^2;- = LL^-P(^^J ' (38) 

and, as we have seen in section 3, there are at least two Hamiltonians associated to Eq. 
which are given by 
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H. 



(1) 



and 



(2) 



2;z 



N2 

= E 

k=l 

E 

k=l 



' 2a 



In 



1 - tanh^ 



« P2zfe 



+ m2gq2zk 



Plzk 

2ms 



exp 



2zk 



m2 



2a 







exp ^ 






m2 y 



(39) 



(40) 



Therefore, one has two different Hamiltonians to describe the same system, H^^^ — Hi.x,y,z + 
H2;x,y + H^l and if^^^ = Hi;x,y,z + H2-x,y + H^], which are written as 



■ 3 2 

EPlij , 
^+mmizj 



,i=l 



2mi 



N2 3 



N2 



+E 

k=l 



P2zk 

2m' 



exp 



1712 



2a 



k=l i=l 



exp 



P2ik 

2m2 

2aq2. 

1712 



exp 



zk 



2aq: 



2ik 



1712 



(41) 



and 
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N 



EPlij , 

j=l lt=l 



^2 3 /n 

j-^ ^ 2m2 V "^2 

k=l 1=1 ^ 



N2 



k=l 



2a 



1 - tanh^ 



+ m2gq2zk 



(42) 

Then, one can calculate for each Hamiltonian the canonical partition function (Toda et al 
1998) which is associated to the same statistical system, 



1 



exp 



1,2 



(43) 



Ni\N2\h^^ 

where [3 is defined as j3 = 1/kT with k being the Boltzman's constant and T being the 
temperature, and the integration is carried out over all the coordinates and linear momenta 
of the two particles. The integration of momenta is carried out in the intervale (— oo, +oo). 
The integration on the transverse coordinates {x,y) is carried out in the intervale [0,L], and 
the integration of the vertical coordinate is carried out in the intervale [0, z] . The partition 
functions for both cases are given by 



NilN2lh^^ \ P J V 



( 



271171'; 



(3mig J \ P 



3N2/2 



e 2a 



2pag 



Erfi 



2a ' •' 



/3gmj 
2a 



N2 



(44) 



and 
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Ni\N2\h^^ V P 



27rmi V^^/^ ^1 _e-/3w\^^ /27rm2\^' /m2\2^2 



(5m2g 



7ra 



V /5 
r 



2a 



n N2 



mlg r + l\ 

\ 2a 2j _ 



The system has two internal energies, U^^^ — —d In Z^^^ / dp, 



and 



5iVi \ 1 



1 _ e-Pmigz 2^ 



(45) 



(46) 



5A^i 



2 

N2m\g 
2a 



+ 2N2 



1 



f/3rr4g\ 
\ 2a J 



Pmlg 



2a 



(47) 



where the function /(/3) (/'(/5) = df{(3)/d(3) has been defined as 



(48) 

and Erfi{x) = —iErf{ix) is the complex error function which can be expressed in the form of 
the Dawson's integral, Erfi{x) = Dawson{x), Dawson{x) — and ip is the 

digamma function, ip{x) — d\nr{x)/dx. Thus, one can have two heat capacity expressions 
for the system, C^^ = dU^'^/dT = -k(3'^dU^y d(3, 

2Ar,J k (1 _ e-™.../3)2 + ^2^/5 [j^ jj^) (49) 
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and 



(2) 



+ 2N2\k 



Nik{migz(3fe-"''<^'f^ 



N2k {m2gzl3f e-'^^s^/s 



+N2k 



2a 



(1- 



2, 



m4g£_ 

2a 



(1) 



2a 



1' " 



(50) 

where ^'^^^ is the trigamma function, ip^^^x) = (P'V{x)/dx^. Figure 1 shows the difference 
\Cy^ — Cy\ as a function oi (3 = 1/kT. As one can see, this difference is not small at 
low temperatures (high /3 values). From /3 lower than about 2100, is higher than Cv^, 
and the situation is reversed for higher values. This difference seems to have an important 
implication related with the ergodic hypothesis (Toda et al 1998). Assuming the validity of 
the hypothesis, one would expect not difference at all on the calculated heat capacities (or 
internal energies) since averaging over the time variable must bring about the same value for 
both Hamiltonians (they represent the same dynamical system). However, averaging over the 
canonical ensemble must be different if the Hamiltonians are different. This ambiguity will 
remain when quantum canonical ensemble is considered (using Eq. ( pTl) and Eq. (1321) ) for 
quantum statical analysis of the system. 



6 Conclusions 

We have shown two constants of motion, two Lagrangians, and two Hamiltonians for a free 
falling particle moving in a media with quadratic velocity dissipative force. These quantities 
describe the same dynamics of the classical system, but their quantization and the associ- 
ated statistical mechanics (for an ensemble of particles) describe two different quantum and 
statistical dynamics. We have showed this at first order in the dissipative parameter and at 
first order in perturbation theory. There is still a point which reamins to to study and has to 
deal with quasi-classical limit. The question is whether or not both quantum Hamiltonians, 
Eq. fl2Tl) . describes the same quasi-classical dynamics {h — > 0) and coincides with the classical 
dynamics in this limit. We will deal with this problem and hope to report some results soon. 
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Figure Captions 

Difference of the heat capacities as a function of /5 = 1/kT for a = 0.01, g = 1, and 
mi/m2 = 0.1 
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